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Abstract 



We examine the long time behaviour of A + B — > reaction diffusion systems with 
initially segregated species A and B. All of our analysis is carried out for arbitrary 
(positive) values of the diffusion constants D^, Dg, and initial concentrations Oq 
and bo of A's and B's. We divide the domain of the partial differential equations 
describing the problem into several regions in which they can be reduced to simpler, 
solvable equations, and we merge the solutions. Thus we derive general formulae for 
the concentration profiles outside the reaction zone, the location of the reaction zone 
center, and the total reaction rate. An asymptotic condition for the reaction front to 
be stationary is also derived. The properties of the reaction layer are studied in the 
mean-field approximation, and we show that not only the scaling exponents, but also 
the scaling functions are independent of Da, Db, a and b . 
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1 Introduction 



The study of the interfacial region formed in diffusion limited A + B — > type reac- 
tions between domains of unlike species has attracted much current interest fl 



A natural way to examine this problem is to prepare a system with the components 
initially segregated along the plane x — 0, and then investigate the spatio-temporal 
evolution of their concentrations pa and p B , and the reaction rate R. Such geome- 
try, first studied by Galfi and Racz f2|, was already investigated by means of various 
methods, including experiments || [|, §], numerical simulations ||, |7], [|, §], analytical 
computations flO], [II], f[2| , scaling [§, [| 0, [bf] and dimensional |~3], [L5| analysis. 



A standard way to treat the initially separated problem analytically is to solve 
the following partial differential equations M 



with the initial state given by 

p A (x,t = 0) 
p B {x,t = 0) 




(2) 



where pA(x,t) and pB(x,t) are the local concentrations of A's and B's, R is the 
reaction rate, H(x) denotes the Heavyside step function, and a , b , Da and D B are 
some positive constants related to the initial concentrations of species A and B and 
their diffusion coefficients respectively. It is customary 0, [T5], 0, [H], |T3|, 



to assume Da = Db = D, which leads to the conclusion that u(x,t) = pa~ Pb obeys 
the readily solvable diffusion equation d t u = Dd\u irrespective of R. Finally some 
form of R must be assumed, and in most cases either the mean field approximation 
R oc PaPb II, Oil HI, mi, or its generalization R oc p™p n B f| ||, [13], |14| was adopted. 



With these assumptions, two fundamental concepts were developed, both referring 
to the long time limit. According to the first one [[|, the long time behavior of the 
system inside the reaction layer can be described with a help of some scaling functions 
Sa, Sb and Sr through 



PA (x,t) oc t^S A (?—^) , (3) 
PB (x,t) oc t^S B (?^M\ , (4) 
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R(x,t) oc t^S R [±-^\ , (5) 

where Xf(t) denotes the point at which the reaction rate R attains its maximal value, 
and exponents a, (3 and 7 are some positive constants given, for R oc p^Pb, by 
7 = l/(m + n + 1), a = | — 7 and /3 = 1 - 7 0. The scaling ansatz is based on 
the assumption that the width w(t) of the reaction layer grows with time as t a with 
a < 1/2, so that in addition to the diffusion length scale Xo ~ VDt, the problem 
possesses also another relevant length scale w oc t a . 

According to the second theory, called the quasistationary approximation [JIB], [TB| , 
the currents J A {t) and Js(i) of particles A and B arriving at the interface layer 
from the two densely occupied domains are changing so slowly, that the relatively 
narrow interface has enough time to equilibrate. To 'equilibrate' means here to reach 
a state completely determined by the current boundary conditions, i. e. by J A and 
Jb- Mathematically this is equivalent to the assumption that the state of the reaction 
zone is entirely given by equations obtained from flU) by replacing their left sides, or 
the time derivatives, with zero. This leads to much simpler equations 



d b %? = R 



(6) 



which are to be solved with the boundary conditions dpA/dx — > —J^it) and ps — > 
as x — ► —00, and pa — > 0, dps/dx — ► Je(t) as 2 — ► +00. The most important 
feature of the quasistationary equations (H) is that they depend only on x, with time 
t being a parameter entering their solutions pA(x,t) and ps{x,t) only through the 
time dependent boundary currents J a and 

It was conjectured by Galfi and Racz |2j that the first of the above assumptions, 
Da = Db, is irrelevant with regard to the long time behavior of the system, the ratio 
Da/D b affecting perhaps the form of the scaling functions S A , Sb and Sr, but not the 
values of exponents a, (5 and 7. This hypothesis was generally accepted after numer- 
ical II and experimental j|] verification. However, there is still no analytical theory 
referring to the general case Da ^ Db- For two reasons this situation arouses some 
anxiety. First, it is practically impossible to find in Nature two species with exactly 
the same diffusion constants. Second, the above mentioned verification encompassed 
only the case where the ratio Da/Db was of order 1, whereas it is known || [17], that 
if one of the diffusion constants is equal zero, the mean-field exponents assume values 
entirely different from those predicted by Galfi and Racz, namely a = 0, (3 = 1/2 and 
7 = 1/4. The aim of this paper is to present such general theory comprising the case 
of any positive diffusion constants Da and D B - 



2 



Unfortunately, we know of only one successful attempt to derive the macroscopic 



form of R from the microscopic properties of the system [11|. Dimensional analysis 



leads to another important conclusion that the mean field approximation should be 



valid only in spaces of dimension higher than d c = 2 [g, [13], [L5| . Therefore our basic 
equation ([!]) might seem useful only for these two sorts of systems for which the form 
of R is known. In our approach, however, we will not need to impose any special 
restriction on the form of R. Instead, we will require that the solutions of (|l|) satisfy 
a few physically justifiable relations. Therefore our theory can be applied even to 
the systems for which the form of R remains unknown, including experiments and 
microscopic models. In such cases verification of our postulates should be far easier 
than the task of finding the exact form of R, let alone solving (H) afterwards. 

The paper is organized as follows. In the next section we will present the assump- 
tions our theory has been founded on, as well as their short physical justification. 
The general theory is formulated in the third section. In the next section we will use 
it to derive and discuss the scaling ansatz in the mean field approximation. The final, 
fifth section is devoted to conclusions. 



2 Assumptions 



We will consider systems which can be described with the Galfi and Racz equations 
(H) and the boundary conditions (§). We will assume that Da, D B , a and b are some 
known positive constants. Our analysis will be based on a few physical assumptions. 

i. At any time t > there exists a unique point Xf(t) at which the reaction term 
R attains its maximal value, and a unique point xo(t) at which DApA{xo,t) — 
D B p B (x ,t) = 0. 

ii. The reaction is concentrated in a region \x — Xf \ ~ w(t) ~ t a with < a < 1/2. 
Outside this region, for x <C Xf — w, there is pa 3> p B , and fori>i/ + rowe 
have pa < Pb- 



iii. The evolution of pa in the region x <C Xf — w can be approximated by 

p A (x,t) = a -C A 



ert\x/y/4D A t) + 1 



(7) 



where Ca is a constant, and erf (x) = 2ir x l 2 Jq exp(— rf) dr/ is the error 
function liT 



Similarly, for x ^> Xf + w, the evolution of p B can be estimated by 
p B (x,t) = b + C B \erf(x/Jw B t) - 1 
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where Cb denotes another constant. Both Ca and Cb depend on the initial 
parameters a®, b , Da and Db- 

iv. The quasistatic approximation is valid in the region — (.D^) 1 ^ 2 "C x <C (Dst) 1 / 2 . 

The first assumption introduces two functions Xf(t) and Xo(t), restricting the 
considerations to the cases where they are uniquely defined. Function Xf identifies 
directly the location of the reaction layer at time t, and Xq is an auxiliary, mathe- 
matical object helpful in examining the behavior of Xf. That xq exists for any t > 
stems from the initial conditions @. As for the second postulate, it was satisfied 
by all the A + B —>■ interfacial systems examined so far. The third assumption 
comes from the observation that, due to postulate ii), in the region x <C x f — w the 
concentration of particles A is expected to be much bigger than that of B's, the latter 
having to cross the whole reaction layer to get there. Therefore, the evolution of 
A's is practically unaffected by B's, and so it should be governed by the standard 
diffusion equation d t pA = Da9 2 Pa- The particular, based on the error function form 
([?]) of its solution was predicted and experimentally confirmed by Koo and Kopelman 
0. Notice also that for any time t such form of pa guarantees that the relation 
lim 3 .^_ 0O p j 4 = a Q implied by the initial conditions (|2|) is also fulfilled. A similar ar- 
gument leads to As for the last postulate, the quasistationary approximation is 
based on the following observation |13|] . The diffusion current of particles arriving at 



the reaction layer is J oc t 1 / 2 , so the characteristic time scale on which this current 
changes is tj ~ (dlog J/dt)^ 1 oc t, whereas the equilibration time of the reaction 
front is Tp ~ w 2 oc t 2a ; therefore a < 1/2 implies that as time goes to infinity, the 
ratio tf/tj goes to 0, validating the quasistatic approximation. 

As we mentioned above, we will not impose any explicit restrictions on the form 
of the macroscopic reaction rate R requiring only that it be consistent with the above 
postulates. However, to investigate the behavior of the A + B — > system inside the 
reaction zone we will need more detailed information about R. Therefore in Section 
4 we will concentrate on the mean- field approximation R oc paPb- 



3 Analysis 

The following observation constitutes the basis of the analysis of our model. For 
sufficiently long time t, at any point x we can employ either assumption iii) or iv) or 
both of them - see Fig. [T| Therefore we can divide the x axis into several regions, and 
in each of them the initial problem of solving (|l]) can be reduced to a much simpler 
one. Then, the overlapping of the domains of applicability of iii) and iv) will enable 
us to merge the solutions. 
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iv 

quasistationary approximation 



in 111 

free diffusion of A's free diffusion of B's 



1 1 1 1 ~1 " 

"°° -VDJ x f -w x f x f + w 00 



Figure 1: Schematic diagram of the regions of applicability of postulates iii) and iv). 
Asymptotically w(t) oc t a < t 1/2 . 

Consider first the region —y/D^t Cx< \fTJ^t. By assumption iv) the system is 
governed here by quasistationary equations (|J). They imply that ^(x, t) = DbPb — 
DaPa satisfies d 2 ^> jdx 1 = 0. Therefore \I/ is linear in x. Let J(t) denote its slope. 
By definition of Xq we have \I/(xo,t) = 0. Thus we arrive at the conclusion that at 
sufficiently long time t, for —y/D^t <i< y/D B t, there is 

DbPb ~ DaPa ~ J(t)(x - x (t)) , (9) 

and so J Ait) = J sit) = J if)- The notation fit) ~ git) means lim^oo fit)/ git) = 1. 

Consider now the region — a/ZJ^ <c x <C Xf — w, so that e = xj — x fulfils 
f < e C t l l 2 . Applying assumption ii) to (|9|) we can approximate the form of pa by 

p A (x,t) « -^^(s-xoCt)) • (10) 

On the other hand, however, by assumption iii), p A can be here as well expressed by 
equation (|7|). So we have 

-DA X J{t){x f it)-x it)-e) (11) 



and 

-L- 

By assumption ii), for any x located outside the reaction layer, the ratio Pa/ Pb 
will either converge to zero, or diverge to infinity as t —>■ oo. However, by definition 
of x , this ratio assumes the constant value D b /Da at x — x . So x must lie inside 
the reaction layer. As its width grows as t a , we conclude that there must exist a 
number 9 such that \xf(t) — xo{t)\ < 6t a . We can see now that in the long time limit 
\xf — xq\ becomes negligibly small compared to e which, in turn, gets negligibly small 



a 



erf 



V^DaI 



erf 



x 



V^DaI, 



+ 1 



DTJit) 



12) 
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compared to t 1 / 2 . Therefore we can drop e on the left hand side of ([11]) and fll2|), and 
Xf — Xq on the r. h. s. of fllTl). After these transformations the asymptotic value of 



the 1. h. s. of (PI) turns out independent of e, whereas the r. h. s. of (|TT]) becomes 
proportional to eJ(t). As e can vary between t a and t 1//2 , we conclude that J(t)e(t) 
goes either to 0, or to oo. The latter case is impossible because ( [TT| ) approximates 
the value of pa which must be finite. In the long time limit we therefore have 



J(t)e(t) 







erff^Ul 
I V^DaI 



and 



J(t)Vt -> C A yjD A /'Kexp^- 



x 



AD A t 



Similar arguments applied to the region ij+wCiC \JD B t lead to 



erf 



and 



J{t)Vt -> C B ^jD B /ir exp - 



AD B t 



It follows from fll4|) and (|T6"|) that in the long time limit 

z/(t)M - C, , 
where Cf is a constant given either by 

C> = 2/D^erf- 1 [K-C A )/C A ] 

or 

C f = 2jD B ~eir 1 i(C B -b )/C B ] . 
Now (|l^), C irri) and ( |i"8|) imply that as time goes to infinity we have 

j(t)Vt - Cj , 

where Cj is another constant given either by 

' C]_\ 



Cj = C A JD A /Tcexp 



(13) 
(14) 

(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



f22^ 
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or 



Cj = C B jD B /7rexp 



AD 



B , 



(23) 



Notice that (E^) is consistent with (|T3|). 

So far we have introduced four constants Ca, Cb, Cf and Cj. The first two of 
them, Ca and Cb, control the asymptotic profile of the majority species outside the 
reaction layer. The third constant, Cf, governs the location of the reaction layer 
center. Finally, through the formula J(t) »s JR(x,t)dx ~ Cj/t 1 ^ 2 , parameter Cj is 
related to the magnitude of the current J(t) of particles entering the reaction layer, 
or, equivalently, the total reaction rate at time t. Due to the form of the initial state 
(0) we expect d x pA < and d x ps > 0, which implies Ca > 0, Cb > and Cj > 0. 



Equations (]TP|), (fZDD, ( |2"2"| ) and (^3|) can be reduced to 



$1 



-C f 



q /DI / C f 



(24) 



where 



$(x) = [1 — erf (x)} exp(x 2 



(25) 



An important feature of $(x) is that it diminishes monotonically from oo to as 
x grows from — oo to oo. This property guarantees that equation (^) always has 
a unique solution Cf = Cf(ao/bo, Da, Db) which, moreover, can be readily found 
numerically. The only problem that can appear while solving (^) numerically is that 
when x is positive, Q(x) is a product of a very small and a very big numbers. For this 
reason, if x is greater than 5, we suggest to use the asymptotic form $(x) ~ l/(^x) 
which comes from the asymptotic properties of the error function erf . 

With Cf computed from (E5), the values of Ca, Cb and Cj can now be calculated 
from fllTf), (|20D and (|22"D. The opposite statement is also true: if we know (e. g. from 
an experiment) the values of Ca, Cb, Cf and Cj, our equations determine uniquely 
the values of ao, bo, Da and Db- 

The immediate consequence of (24) is that the sign of Cf is determined by the 
sign of qq^Da/ (bo^Ds) — 1. In particular we conclude that 



C f = ^ a sjD A = boy Db . (26) 

This formula is important for planning experiments, as it clarifies the way the initial 
concentrations of the species should be chosen in order to have the reaction layer 
move asymptotically as slowly as possible. Condition (|26|) is consistent with that of 
Jiang and Ebner's |J who, by numerical examination of the mean field approximation 
R oc PaPb, found a stronger relation Xf = <^=^> clo^/Da = bo^Ds- Our general 
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formula, derived for any reaction term R, implies only that with this particular choice 
of the initial parameters the function Xf cannot be changing as fast as t 1//2 . An 
example of a system where Cf — and Xf(t) oc t a was investigated in \TA . 

Equation (^4j) enables us also to observe a striking similarity between the long 



and short time behavior of Xf. According to [K|, in the short time limit the reaction 
term does not affect the solutions of (|l|), and so pa and pB assume the same forms as 
in the readily solvable case R = 0. The point Xf can be then found as the point at 
which dR/dx = 0. For R oc p™p B such procedure yields lim^ Xjj\ft = Co, where 
Co can be found from the relation very similar to that of (El) 



$1 



Cn 



$1 



-C 



(27) 



4 The reaction layer 

In the previous section we carried out our analysis without imposing any restrictions 
on the form of the macroscopic reaction term R. As we now proceed to examine 
the asymptotic properties of the reaction layer, we will obviously need more specific 
information about R. Therefore we will concentrate on the mean field approximation 
R = kpAPB, k = const., still allowing ao, bo, Da and Db to take any positive values. 

By assumption iv) we expect that in the region -{D A t) 1/2 < x < {Dst) 1 ' 2 we can 
apply the quasistatic approximation equations (Q). Let pa(x, t) and pB(x,t) denote 
their solutions for some values of Da, Db, Xo(t) and J{t). By the following linear 
transformation we introduce two new functions of a single variable Pa{%) and p~b{%) 

p A {x,t) = r) A (t)p A [(x-Xo(t))/w(t)} ) 

\ (28) 
p B {x,t) = r)B{t)p B [(x - x (t))/w(t)) J 

where 

/ D \ 1 / 3 / C \ 2//3 
VA (t) = J(t)w(t)/D A = (-f) (-^j r 1/3 , (30) 

/ n \ x /3 / n x 2/3 
VB (t) - J(t)w(t)/D B = (^) (g£) . (31) 

Denoting R(x) = pa(x)p~b(x) we have also 

R{x,t) = kp A pB = Cj^iDADB^^k^H-^RUx - x )/w(t)] . (32) 
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The essential property of p~a{x) and /5b (x) is that they constitute the particular 
solution to equations (H) with D A = D B = J = k = 1 and x = 0. Therefore, by sym- 
metry, R(x) assumes its maximal value for x = 0, and so equation fl32|) implies that R 
attains the maximal value at x = xq. In the long time limit we can therefore identify 
Xf with xq. Comparing now fl28|) and fl32|) with the scaling ansatz (|3|) - (|5|) we see 
that we can also identify Pa with Sa, Pb with S B , and R with Therefore not only 
are the above formulae consistent with Galfi and Racz's scaling ansatz, but through 
Cj(a ,b , Da, D B ) they also exactly relate the quantities of physical importance (e. 
g. w(t)) to the parameters of the system (Da, Db, cl , 6 and k). 

Because (|^) can be applied to systems with any positive values of 'external' 
parameters a , b , Da, Db and k, we arrive at the conclusion that the long time 
evolution of initially segregated A + B — > systems is even more universal that it 
was predicted by Galfi and Racz; namely, not only the scaling exponents, but also the 
form of the scaling functions does not depend on the external parameters. Therefore, 
to find the scaling properties of the reaction layer it is sufficient to concentrate on the 
simplest, symmetric case Da = Db and = bo- 

Notice that we have achieved these results by means of a simple, linear transfor- 
mation In this way we took advantage of the very feature of equations (]!]) and 
(P) that prevents them from being solved analytically - nonlinearity. 

The above analysis is straightforward and can be easily generalized for many other 
reaction terms R. In particular, for R = kp^p B , with k = const and m, n being any 
(positive) real numbers, the following relation should be used instead of ( |29|) 



w m+n+L = D™D n B k- L J L - m ~ n . (33) 

This formula, together with (|2TD, ( PU[ ) and Q3"T|), generalizes the scaling theory of 
Cornell et al [[| for the case of any positive ao, b , Da, Db, m and n. 



5 Conclusions 

We have investigated the long time behavior of the concentrations pa and ps of phases 
A and B in the Galfi and Racz's problem. Our analysis is the first analytical attempt 
to consider it in the general case of arbitrary positive initial concentrations a$ and b , 
and diffusion constants Da and Db of A's and B's. 

Our approach is very general, as it does not impose any restrictions on the form 
of the macroscopic reaction rate R. Instead, it is based on the assumption that in the 
long time limit pA and ps satisfy a few physically justifiable relations. Therefore our 
theory can be applied to various systems, including those for which the form of the 
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macroscopic reaction rate R remains unknown. Another peculiar feature of our theory 
is that, unlike most of previous studies, it does not concentrate on the investigation 
of the reaction layer only, but takes into account the properties of the whole, infinite 
system. 

In this way we managed to derive general formulae for the concentration profiles of 
the majority species outside the reaction layer, the location of the layer, and the total 
reaction rate. It is interesting to notice that these quantities turned out independent 
of R. We also derived analytically Jiang and Ebner's condition for the reaction front to 
be asymptotically stationary. This relation also turned out independent of R. These 



results correspond to the recent findings based on dimensional analysis |9], [13], |15 
according to which the scaling properties of the reaction layer are independent of the 
form of R. 

Next we derived the general scaling ansatz for the mean field approximation. We 
gave the formulae which exactly relate some quantities of physical importance, (e. g. 
the width w of the reaction layer) to the external parameters of the system a , b , 
Db and k. It turned out that not only the scaling exponents, but also the forms of 
the scaling functions are independent of the values of these parameters. This justifies 
the customary approach of examining the properties of the reaction layer only in the 
simplest, symmetric case ao = b and Da = Db- 

Our work suggests also that the behavior of the reaction-diffusion system can 
be understood as a subtle interplay between two scaling regimes. The first one 
is valid far from the reaction zone, where the densities of particles A and B as- 
sume the scaling forms typical of purely diffusive systems: pa{x, t) ~ ^(x/t 1 / 2 ) and 
Ps(x,t) ~ b{x /t 1 / 2 ). These scaling laws determine also the location of the point 
Xf(t) of the maximal reaction, and the magnitude of the current J(t) of the particles 
entering the reaction zone. However, at Xf the spatial derivatives of and ^b suffer 
discontinuity. Therefore in the vicinity of Xf a new form of scaling develops, and pa 
and pb assume the form pA(x,t) = SA(x/t a ) and = S^x/t") with a < 1/2. 

Although we confined our considerations only to the long time limit, it would be 
interesting to combine our results with those of Taitelbaum et al [1(J for the short and 



intermediate times. We believe that the striking similarity between equations (|24| ) 
and (|27|) is not accidental and should lead to a general theory comprising the short, 
intermediate and long time limit. The first attempt in this direction has already been 



made fl8fl 



Notice also that the quasistationary approximation leads to new definitions of 
'short', 'intermediate' and 'long' time regimes. Namely, we can define them as the 
time intervals in which the reaction term, in the vicinity of Xf, is vanishingly small 
compared to the time derivative ('short time'); or the interval in which they are 
of similar magnitude ('intermediate time'); or the interval in which it is the time 



10 



derivative that can be neglected ('long time'). 

Another interesting problem concerns the limit Da — > with other external pa- 
rameters fixed. In this limit the scaling exponents (in the mean-field approximation) 
are expected to change from a = 1/6, f3 = 2/3 to a — 0, /3 — 1/2. The paper in 
which this problem is examined within the framework of the presented here theory 
is under preparation. We will mention here only that as Da goes to 0, the time at 
which the system reaches the long time regime goes to infinity so the case Da = 
can be considered as the case where the system always remains in the 'intermediate' 
time regime. 
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